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O ! We briefly review the existing \jf(2S) data taken at RHIC, the Tevatron and the LHC. We system- 
ic | atically compare them with colour-singlet-model predictions as a function of the center-of-mass 
energy, of the quarkonium rapidity and of the quarkonium transverse momentum. The overall 
i— — i agreement is good except for large transverse momenta. This points at the existence of large 
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NNLO corrections or points at colour-octet dominance. 
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1. Introduction 

The physics of quarkonium production in high-energy hadron collisions has triggered much 
investigation and debate since the preliminary data of prompt y(2S) production at the Tevatron 
were released by the CDF collaboration in 1994 [1]. It uncovered an obvious discrepancy between 
the measured Pj spectrum and the one predicted by the simple application of pQCD, dubbed the 
colour-singlet model (CSM): the spectrum of the data was harder with a clear excess of the y{2S) 
at large Pj. This observation was confirmed in the final publication in 1997 [2] and similar obser- 
vations were also made for the J/y -once the feed-down from the % c could be subtracted [3]- and 
for the T family [4] -albeit to a less extent. 

Until not loo long ago, the most popular explanation of these discrepancies was provided 
by an enhancement of the yield due to non-perturbative transitions between the colour-octet state 
produced at short distances and the colour-singlet mesonic states (for reviews see [5]). However, it 
is clear nowadays that one cannot neglect the 0$ and aj corrections to the CSM [6] if one wants 
to confront the prediction of the CSM with the data for the Pj spectrum of J/y and T produced in 
high-energy hadron collisions [7, 8, 9, 10, 11]. In fact, this indicates that a factorised description 
of high-Pr quarkonia beyond leading power [12] may be more suitable than NRQCD [13]. 

QCD corrections also have a dramatic impact on the predictions of polarisation observables. 
At LO, the yield of inclusive J/y and T is transverse -in the helicity frame-, whereas it "becomes" 
longitudinal at NLO at large Pj [9, 10]. This is also true for instance when they are produced in 
association with a photon [14, 15]. The sole known case where the quarkonium polarisation is not 
altered by NLO corrections is when they are produced with a Z-boson [16]. 

On the other hand, in recent works [17, 18], we have shown that the CSM alone is sufficient to 
account for the magnitude of the Pj -integrated cross section and its dependence in rapidity, da/dy, 
at RHIC, Tevatron and LHC energies. In other words, there is no need for additional contributions 
at low Pj and the energy dependence is well reproduced. 

We consider here the \fs, Pj and y dependences of the y(2S) yield at collider energies. We 
present various comparisons between the existing data and the yield at LO, NLO and sometimes 
including some dominant contributions at aj (NNLO*) when possible. It should be stressed that the 
NLO code for the CSM at high energies and low Pj is not stable likely owing to large contributions 
of the loop corrections which can become negative at low Pj. NLO and NNLO* are thus only 
computed for sufficiently large Pj. 

After having briefly described how the y(2S) cross sections are evaluated in the CSM, we 
compare them with the existing data as a function of the cms energy, of the quarkonium rapidity 
and of its transverse momentum. As we shall see, the overall agreement is good except for large 
transverse momenta. 

2. Cross-section evaluation in the CSM 

Details about the evaluation of quarkonium production cross section at LO in the CSM can 
be found in [6]. As regards the evaluation of the yield at NLO accuracy, we have used the code 
of Campbell, Maltoni and Tramontano [7]. In the CSM, the cross section to produce a y(2S) is 
obtained along the same lines as for the J/y ground state with the mere change of the value of 



2 



y/(2S) production in proton-proton collisions at RHIC, Tevatron and LHC energies 



J. P. Lansberg 



the wave function at the origin R(0), which is meant to account for all the non-perturbative and 
relativistic effects. For the J/yf, \R{§)\ 2 can be taken as 1.01 GeV 3 ; here, we took 0.67 GeV 3 for 
the y{2S) [17]. 

The expected impact of NNLO QCD corrections for increasing Pj is investigated by evaluating 
what we call the NNLO* yield. A complete discussion can be found in [10, 11]. In a few words, 
we anticipate that the NNLO* yield encompasses the kinematically-enhanced topologies which 
open up at aj with a Pj fall off of do jdPj. Let us emphasise that we do not foresee significant 
modifications of the Pj dependence at N 3 LO and further. At NNLO, the Pj A fall-off of the new 
NNLO topologies is the slowest possible. Above aj, the common wisdom about the decreasing 
impact of further QCD corrections would then hold. One expects a K factor multiplying the NNLO 
yield to be independent of Pj and to be of the order of unity. It would be worrisome to find 
out a further enhancement between the NNLO and N 3 LO results by an order of magnitude. As 
in [10, 1 1], it is evaluated thanks to a slightly tuned version of the automated code MADONIA [19]. 

The uncertainty bands at LO and NLO are obtained from the combined variations of the charm- 
quark mass (m c = 1 .5 ± 0. 1 GeV), the factorisation ±ip and the renormalisation [Ir scales chosen in 
the couples ((0.75,0.75); (1, 1); (1,2); (2, 1); (2,2)) x m T with m\ =4m 2 Q + P%. The band for the 
NNLO* is obtained using a combined variation of m c , 0.5mr < Hr = Hf < 2mr and 2.25 < s™ m < 

9.00 GeV 2 . We have used the LO set CTEQ6_L and the NLO set CTEQ6_M [20] and have taken 
l^v(25)(°)l 2 = °- 67 Gey3 and Br(v/-(2S) -> £+£-) = 0.0075. 

3. Result and comparison with existing data 

3.1 Pj -integrated yields 

In Fig. 1, we compare the y/s dependence of do/dy\ y= o obtained at LO in the CSM and that 
obtained from the PHENIX data (y/s = 200 GeV) [21], the CDF data 1 (y/s = 1960 GeV) [22] and 
rescaled LHCb data (y/s = 7000 GeV) [23]. As one has obtained for the J/y and the T(nS) [18], 
one does not observe any surplus w.r.t. to the CSM predictions. In fact, at high energy, the LO 
yield tends to be above the experimental data. 

As mentioned above, we are unfortunately not able to give, for the time being, corresponding 
predictions at NLO since these are not well behaved at high energies. This is maybe due to the 
contributions of the loop corrections which change sign for Pj of the order of the quark masses. 
This calls for the resummation of initial state radiations (ISR). Yet, at RHIC energies, we have 
observed in [17] that the NLO yield lies in the lower range of the LO uncertainties (see Fig. 2a). It 
is therefore sound to expect that the slight overestimation of the LO w.r.t. the present data would 
be reduced at NLO once ISR can be resummed and stable results can be presented. Fig. 2a and 
Fig. 2b shows the rapidity dependence at Tevatron and LHC energies. Note that, for now, there is 
no measurement of the y(2S) yield for low enough Pt to derive a Pj -integrated cross section for 
y ~ 0. 

'To be precise, the CDF data only covers the region Pj > 2 GeV, for which (j(\y\ < 0.6, Pj > 2 GeV) x Br = 
2.6 ± 0. 1 nb. We have assumed the same Pj dependence as the inclusive J / \f/ for Pj > 2 GeV [25] and we have obtained 
<j(|y| < 0.6, Pj < 2 GeV) x Br = 2.1 ±0.1 nb. On the way, it is instructive to keep in mind that 45% of the J/y yield 
lies below Pj = 2 GeV. 
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Figure 1: Data vs LO CSM for the Pj -integrated y{2S) production cross section at y = as a 
function of the cms energy. The forward LHCb data are rescaled assuming the same y dependence 
as the LO CSM, i.e. using a factor 1.5 for do / dy\ y= o / (do / dy) |2.o<y<4.5- Data are from [21, 22, 23]. 




(a) RHIC (b) LHCb 

Figure 2: Data vs (N)LO CSM for the -integrated direct-quarkonium production cross section as 
a function of the y(2S) rapidity. Data are from [21, 23]. 

3.2 P T -differential yields 

Now, we move onto the discussion of the Pj dependence of the cross section. As we discussed 
in the introduction, it is clear that the LO CSM cannot be sufficient since the leading (and even 
sub-leading) Pj topologies are missing. Except for low Pj, where the results are not stable, the 
NLO yield can nowadays be computed easily and compared to the data. It is shown at RHIC and 
Tevatron energies in Fig. 3a and Fig. 3b; one clearly sees a different Pj dependence w.r.t. the LO 
yield. It is clearly harder and the discrepancy with the experimental data is reduced, altough it is 
still significant when Pj gets large. 

Indeed, at large Pj, one expects the Pj kinematical enhancement of some NNLO topologies 
to take over the extra a s suppression w.r.t. the NLO. This is why the NNLO* bands in Fig. 3a- 
3d show a Pj spectrum which is even harder and much closer to the data. If the CSM at NNLO 
is indeed the physically relevant production mechanism at large Pj, namely that the charmonia 
are produced at large Pj along with two hard partons (in a sense, two jets), the corresponding 
theoretical predictions would be involving. Even with a full NNLO, one would have to face large 
theoretical uncertainties due to the factorisation scale through five powers of CC S . In the coming 
years, it is not clear that one could make precise and definitive comparisons between data and 
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theory as far as the Pj dependence of the yield is concerned. It may thus be more fruitful to also 
analyse new observables such as associated production. 
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Figure 3: y(2S) data vs CSM predictions at various orders for the Pj dependence of the differential 
cross section at RHIC, Tevatron and LHC energies. Data are from [21, 22, 23, 24]. 



For the time being, the most one can tell from the data-theory comparison in the framework 
of the CSM and the approximate NNLO* is that below, say Pj = 15 GeV, the upper uncertainty 
band of the NNLO* agrees with the data - the same observation also holds for the J/yr [1 1, 26, 27] 
while, for the T(nS), the agreement is clearly much better [10, 28]. Whether this has some physical 
meaning is a question that would only be answered once one has a full NNLO computation and/or 
once the LHC experimental collaborations release polarisation measurement for prompt y(2S) 
which would be precise enough to rule out the predictions from either the CSM or from the CO 
dominance. 



4. Conclusion 

We have compared the existing data of y(2S) production at RHIC, the Tevatron and the LHC 
with predictions from the CSM for the Pj, sfs, y dependences of the yield. The %fs and y de- 
pendences are well reproduced at LO and NLO, where it is available. It is therefore worth re- 
investigating [29, 30, 31] the possibilities to constrain gluon PDFs with low Pj quarkonium data. 

As regards the Pj differential cross section, the upper bound of the NNLO* CSM predictions 
of the Pj differential cross section is very close to the experimental data, as previously found for 
T [10, 28] and for J/y [11, 26, 27]. However, the NNLO* evaluation is not a complete NNLO 
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calculation. It is affected by logarithms of an infrared cut-off and its effect might not vanish with 
increasing Pj as quickly as one has anticipated. A full NNLO evaluation of the cross section in the 
CSM is therefore awaited for. 
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